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1^ , Abstract. I review recent work ([TJ [2]) on holographic superconductivity with Einstein-Gauss- 

■ Bonnet gravity, and show how the critical temperature of the superconductor depends on 

both gravitational backreaction and the Gauss-Bonnet parameter, using both analytic and 
numerical arguments. I also review computations of the conductivity, finding the energy gap, 
and demonstrating that there is no universal gap ratio, ujg/Tc, for these superconductors. 
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, ■ 1. Introduction 

Holographic superconductivity (see O SI [5] for reviews) is a fascinating idea which proposes to 
describe strongly coupled high Tc superconductors via a classical gravitational system using the 
gauge gravity correspondence, [6]. The bulk theory typically has gauge and a charged scalar 
V/^ \ field, with a black hole providing a finite temperature. Choosing an appropriate scalar potential, 

•/^ ' such as a negative mass-squared above the Breitenlohner- Friedman bound, [7J , allows the scalar 

to condense out of the vacuum near the black hole horizon, if the curvature is sufficiently large 
CnI ; there |8|. This occurs for low temperature black holes, hence the analogy to superconductivity. 

The superconducting phase corresponds to a "hairy" black hole, where the condensation of 
the scalar field out of its symmetric state screens the charge of the black hole. Typical no-hair 
theorems [9] do not apply in this case as the scalar potential is not positive definite. Near 
. the boundary, the scalar has a power law fall-off, and the coefficient of this fall-off can be 

r> \ interpreted as a condensate in the boundary theory. Following the behaviour of this condensate 

■ as the temperature drops reveals the typical behaviour for an order parameter governing the 

superconducting phase transition, and has been observed in a wide variety of bulk thoeries, 
including varying the scalar mass and potential, varying the number of spacetime dimensions, 
as well as having magnetic fields present and the stability of the system [HI [lOl [TTl [121 ttH [13 [E] • 
Similar models have also been found from embeddings in string theory, [T7], and it is therefore 
of interest to consider more general stringy aspects of these models. 

In this presentation, I review the results of [1] and [2], which explored adding the higher 
curvature Gauss-Bonnet (GB) term, [IE], to the gravitational action. I will also update on 
the talk where appropriate, and add a few new results and comments. In brief, many of 
the qualitative features of the holographic superconductors are stable under higher curvature 
corrections, however, there are quantitative differences, in particular, a universality relation for 
an energy gap in the superconductor was found not to be substantiated even in the Einstein 
limit. (See also [19j for other work on GB holographic superconductors.) 

I first review the bulk superconductor, deriving some general analytical results and remarking 
on the probe superconductor. Then I summarize the numerical results, showing how the 



numerical results follow closely an analytical bound, revealing an interesting subtlety in the 
response of the superconductor to the Gauss-Bonnet parameter a. Finally, I review the results 
of [1] and [2] on conductivity. 



2. The bulk superconductor 

The general action is that of an Einstein-Gauss-Bonnet (EGB) gravitational action coupled to 
a massive charged complex scalar field and a U(l) gauge field: 
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where g is the determinant of the metric, and Rabcd, Rab and R are the Riemann curvature 
tensor, Ricci tensor, and the Ricci scalar, respectively. The Gauss-Bonnet coupling constant is 
taken in the range a E [0,L2/4], and the cosmological constant is given in terms of a length 
scale, L, which for Einstein gravity (a = 0) corresponds to the adS curvature. The Planck scale 
is set by = SirG^, q is the charge, and the (squared) mass, of the scalar field. 
The equations of motion satisfied by the bulk system are: 
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Tab is the matter energy momentum tensor 
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and "Da = — iqAa is the gauge covariant derivative. 

To examine holographic superconductivity, we allow for a bulk black hole in order to have a 
system at finite temperature: 
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is the black hole gravitational potential in the absence of backreaction, and 
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is the actual curvature of the adS spacetime, which is renormalized away from the cosmological 
constant scale, L, once a is nonzero. The temperature of the background is: 
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which is equal to rj^/irL? in the absence of backreaction. 

Taking a static ansatz, = (0(r),O...) and ip = tp{r) (which can be taken to be real), the 
equations of motion become 
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where a prime denotes derivative with respect to r. These equations have several scaling 
symmetries, as noted in [2j, which are used to set L = Q = q = 1. Note that by fixing 
Q = 1, the charge parameter is kept fixed in all numerical computations. 

(i) r — ar, t, — >• at, ax^, L — )• aL, q ^ q/a, a ^ a^a, A — )• aA. 

(ii) r ^ br, t ^ t/b, — )- x^b, f — )- 6^/, -)• 

(iii) — )• C(/), — ^ c^i 9 — ^ ^/c, — ?• l(? . 

The horizon is defined in general by f{r-^-) = 0, and demanding regularity of the solution at 
both the horizon and boundary gives the following boundary conditions: 
• Horizon: 
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where A± = 2 ± ^JA + rn?IJl for a general mass m?. We choose C_ = 0, then P and C+ are 
fixed by consistency with the near horizon solution. According to the adS/CFT correspondence, 
we can interpret (Oa^) = C+, where Oa,,. is the operator with the conformal dimension A+ 
dual to the scalar field. 

In [1] we chose to set the mass of the scalar field as = — 3/L^, so that the mass remained 
the same as a was varied. However, the variation of the effective asymptotic AdS curvature, 
d?]), with a relative to L means that this mass actually increases (i.e. becomes less negative) 
with respect to the asymptotic AdS scale, and therefore the dimension of the boundary operator 
corresponding to C+ varies with a. Fixing = — 3/Lgg relative to the asymptotic AdS 
scale, avoids this problem, however the mass now varies with respect to the physical mass and 
temperature of the black hole as a varies. See [20j for a detailed examination of the effect of 
scalar mass on the holographic superconductor. 



2.1. Analytic hounds 

While the equations of motion ([9l ll2p require a numerical solution in general, even in the absence 
of backreaction, near the critical temperature we can glean a measure of analytic information 
using the (uncondensed) charged black hole solution, [21j : 
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and linearizing in the scalar condensate. In the above, Q is the charge of the black hole (up to 
a geometrical factor of 47r), and r_|_ is the event horizon, which determines the "ADM" mass of 
the black hole |22j . To avoid a naked singularity, we restrict the parameter range as a < 

Linearizing in the scalar condensate leaves the background solution for the metric and gauge 
field unchanged to leading order, and we can focus on pO|) with / and (p taking their background 
values. Analytic bounds on the critical temperature can be obtained by looking for simple 
relations that must hold for the existence of a nontrivial scalar solution. First consider the 
variable Xn = r'^ip, which satisfies to leading order: 
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To get an upper bound, let n = 2, and consider the general properties of a solution. At the 
horizon, 
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which is positive for small k^Q'^ (taking X2(r+) > without loss of generality). Since 
X2 ~ l/r2+^ 
that 
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as r — )• oo, the solution must have a maximum for some r, which requires 
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at this point. An examination of when this combination is never positive provides an upper 
bound for Tc- (A different n was used in yy, which sufficed for the probe limit although it led 
to a looser bound. For consistency and comparison between the scalar masses, this new bound 
is used here.) 

We can also obtain a lower bound by considering n = 3. Manipulating ()2ip shows that if a 
solution exists, then the integral 
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is negative. Note that negativity of this integral does not imply existence of a solution to the 
linearized equation near Tc, it is simply a necessary condition. Since this integral is always 
negative at large T, and positive as T — >• (for < 0.4), observing where it changes sign 
provides a lower bound on T^.. This bound was found to give an extremely reliable indicator of 
Tc as computed numerically. 



2.2. Analytic expansion 

In instances where the exact solution is only known numerically, it can be useful to derive 
an analytic approximation to the exact solution using a technique of matching. The solution is 
Taylor expanded around the horizon and boundary, then matched at a midpoint. This technique 
was used in [1] to get an approximation to Tc, although the plot in [Ij was inaccurate. In brief, 
we first transform the radial variable to z = Vj^/r^ and expand the near horizon solution of a 
generic field Y as: 

Y c^YQ-Yi{l-z) + \Y2{l-zf (25) 

The equations of motion give: 
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where /i, /2, and vi are given by (|llll2p . and depend on </)o and -00 • Since the asymptotic 
solutions for the gauge and scalar field given by ()17p also have two unknowns, matching these 
fields and their derivatives at some intermediate point zq gives four equations which determine 
these four unknowns: 



P-^ = -4>i{l-z,) + ^{l-z,f (29) 



-2Qzo 



'^2 



(1 - zo) (30) 



= ^Po-Ml-zo) + t^{l-zo? (31) 



2 



A+^^±4r- = i^i-Mi-zo) (32) 

r+ + 

Clearly, given the complexity of the expressions (j26|) -(|28|). which themselves contain fi, this is a 
rather involved algebraic process in the general backreacting and non-critical case. Nonetheless, 
the system is tractable analytically if we simply focus on deriving the critical temperature. In 
this case, we once again can take the background forms for <j), f and z^, which considerably 
reduce the complexity of the relations: (p9]) and ([30]) are automatically satisfied, = 0, ([12]) 
gives 

/i = "4+^ 

/2 = + ^ + (34) 

and consistency of the ip equations (f3T]) and ([32]) implies 

2A+Vo - 2 (A+ + (1 - A+)zo) V'l + (A+ - 2(A+ - l)zo + (A+ - 2)zi) = (35) 

Although this is a cumbersome expression for nonzero k^, it is just a quadratic for (pf/r'^ = 
AQ'^ /r^, and a linear relation in the probe limit, which in either case can be straightforwardly 




Figure 1. A plot of the critical temperature as a function of a for the scalar mass m? = —"ijl?. 
(a) shows a detailed plot of the probe limit, with upper and lower bounds shown as solid black 
lines, the matching method value as a dotted line, and numerical data points, (b) shows a plot 
of the upper and lower bounds together with exact numerical results for the probe limit in black 
(dotted lines and diamond data points), k} = 0.05 in blue (dashed lines and square data points), 
and = 0.1 in red (solid lines and circular data points). 



solved and the critical temperature obtained from 
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It only remains to choose a value for zq. In [T] the value zq = 1/2 was chosen to be specific, 
however, this choice does not take into account the variation of the effective adS length scale 
with a. In essence, the matching point should be chosen to represent where the asymptotic 
expansion / ~ r^/z^Lg turns over into the near horizon expansion. This obviously depends 
on the relative ratio of r+ to Lg, although as a — )• this point moves rather significantly 

towards the boundary. A good approximation to the departure from the asymptotic regime is 
Zq ~ a/L^ + 1/2, which is the value used in figure [T]||a) . It must be noted however, that there is 
a degree of arbitrariness in this matching process. Although the choice of zq attempts to take 
into account the likely variation of the domains of validity of each expansion with a, this is a 
rather ad hoc process, and as backreaction is switched on, the matching procedure becomes less 
and less reliable. It is included here for completeness, however by far the best analytic guide to 
critical temperature is the lower bound. 



2.3. Discussion of analytic results 

By examining the behaviour of a putative condensed scalar in the bulk, it is possible to derive 
upper and lower bounds to the critical temperature. A matching method can also be used to 
get an approximate solution for the scalar, however, this method (as well as the upper bound) 
becomes unreliable as the backreaction is increased. On the other hand, the lower bound appears 
to be extremely representative of the behaviour of the actual critical temperature, and works 
well for a significant range of backreaction. Figures [1] and [2] show the analytical results, and will 
be discussed further together with numerical data in the next section. 



3. Backreacting superconductors 

In order to find the actual behaviour of the bulk superconductor, in [1] and we integrated 
([9l-ll2p numerically. As already stated, we took L = Q = q = 1, and varied r+ to study how 
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Figure 2. A plot of the critical temperature as a function of a for a selection of k? . The 
analytic bounds are shown as lines and the numerical data as points. Respectively: = is 
shown in black, with solid lines and triangular data points; = 0.05 has blue dashed lines and 
diamonds; k? = 0.1 has a grey dotted line and squares; = 0.2 has a red dot-dash line with 
circular data points. The lower bound is shown for all values, but the upper bound is shown 
only for the lowest two values of k^, as they overlap significantly with the other data and confuse 
the plot. 



the system reacted to varying temperature. From these numerical solutions, we can deduce the 
critical temperature, determine the behaviour of the condensate with temperature, and compute 
the conductivity of the holographic superconductor as discussed in subsection 13.31 

3.1. The dependence ofTc on a and . 

The analytic bounds for Tc have been plotted in figures [1] and [2] together with the exact values of 
Tc obtained by numerical computation for both masses of the scalar field. For = — 3/L^, both 
bounds are shown, as well as the matching method result for the probe limit. For = — 3/Lg, 
the upper bound has only been shown for < 0.05, as above this value it becomes less predictive 
and clutters the plot, and indeed beyond ~ 0.2 (corresponding to (7 ~ 2.25 in the notation of 
[101) it ceases to have quantitative value for any a. The lower bound on the other hand becomes 
successively more accurate as the values of are stepped up, and gives a very good quantitative 
guide to the behaviour of Tc as we vary a and k^. 

It is easy to see that in all cases, the effect of backreaction is to decrease Tc and thus make 
condensation harder. We can see this by using Gubser's rough argument, [8], that the effective 
scalar mass m^g = tv? — (f' (t>{r)'^ / f {r) becomes more negative as backreaction is turned up. 
Essentially, the effect of backreaction is that the condensation of the scalar field not only screens 
the charge of the black hole, but also its mass, as the scalar and gauge fields now contribute to 
the ADM mass. This means that for a given charge and temperature, the radius of the black 
hole is increased, which makes it harder for the scalar to condense. 

One very interesting feature clearly exhibited in the bounds is the turning point in Tc as 
a function of a for the mass m? = —3/Lg. In the probe case, this occurs very near the 
Chern-Simons limit a = and is barely perceptible in the numerical data, however once 

backreaction is switched on, the minimum becomes much more pronounced, and indeed for large 
backreaction (k^ = 0.2) the Chern-Simons limit is showing a considerable enhancement of Tc 
over the typical values for lower a. 




Figure 3. Two plots of the condensate as a function of temperature for = —Zjl? (left) and 
rr? = — 3/Lg (right). The plot on the left shows several values of a in the probe limit: a = 
in red, the lowest line, then a = 0.1 in brown, 0.2 in blue and 0.25 in green. The plot on the 
right shows = in solid lines, and = 0.1 as dotted lines. The black plot is a = 0, green is 
a = 0.125 and magenta is a = 0.25. 



3.2. The scalar condensate. 

The numerical results allow a full exploration of the dependence of the condensate on 
temperature, and we expect to see the condensate switching on sharply at the critical 
temperature, then saturating at low temperatures, and this is indeed what is seen. Figure 
[3] shows (0/^^)^^'^+ as a function of temperature for both scalar masses. 

Each line in the plot shows the characteristic curve of the condensate (Oa+)^^'^^ condensing 
at some critical temperature. For = — the curves are shown in the probe limit. While 
the critical temperature does not display the minimum in a, it can nonetheless be seen that the 
Chern-Simons limit does look rather different, and the scalar field saturates far more slowly as 
the temperature drops. For m? = — S/Lg, three different values of a (0,0.125,0.25) and two 
different values of (0,0.1) are chosen to display the features of the system. Both plots are 
shown with the curves normalized by Tg. In the latter plot, the effect of is to increase the 
height of these graphs, in spite of the fact that the raw data tends to have lower values of {O3). 
This is clearly because the most significant impact of increasing gravitational backreaction is 
that the critical temperature of the system is lowered. In this case the effect of backreaction is 
extremely marked, with the condensate varying more widely with a. 



3. 3. Conductivity 

In |llj . Horowitz and Roberts observed an interesting phenomenon for the conductivity of the 
boundary theory. They considered the model ([1]) in the probe Einstein limit for a range of 
different bulk scalar masses, and on computing the conductivity found an apparent universal 
relation 
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with deviations of less than 8 %. 

Conductivity is conventionally expressed as the current density response to an applied electric 
field: 

a = |. (38) 



The boundary four-current is dual to the bulk field A^, thus we must consider perturbations 
of A^ to compute the conductivity. The full perturbation problem must take into account 
perturbations of the metric [hu and hri) and the perturbation of the Gauss-Bonnet tensor 
Hti, Hri- After some algebra, these reduce to: 

2- •• L^e^" f a{2v'f + f)\ _ 2K^r'^Ai(l)' , , 

^[r/e^^^]'-y^ + ^AA.-|,VA. + -^(/.;.-\^ = (40) 

where h^^ is the perturbation of the metric tensor, and Ai is the perturbation of the gauge field, 
which has only spatial components. Writing j4j(t,r, x*) = A[r)e^^'^~^'^''ei, and setting k = 0, 
p9]) can be integrated and substituted in (jlO|) to obtain: 
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This is solved under the physically imposed boundary condition of no outgoing radiation at 
the horizon: 

A{r) ~ fir)-'^ , (42) 

where T+ is the temperature. In the asymptotic adS region (r — t- oo), the general solution takes 
the form 

. 02 aoLju}"^ ^ r 

^ = ^o + ^ + ^^logj (43) 

where oq and 02 are integration constants. Note there is an arbitrariness of scale in the 
logarithmic term, as pointed out in [11], however, this is related to an arbitrariness in the 
holographic renormalization process; see Appendix A of [2] for a full computation of the 
conductivity and discussion of this renormalization scale. After careful computation using the 
method of Skenderis [23], the conductivity of the EGB system is found to be [2]: 

Note that the imaginary term linear in a; has an arbitrariness of scale from the counterterm 
subtraction, and this is capitalized in the presentation of numerical data, where a suitable 
renormalization scale is chosen to make the features of the plot clearest. 

Figure U] shows the real and imaginary parts of the conductivity as a function of uj/Tc, 
calculated for the scalar mass m? = — 3/L^ in the probe limit and for a selection of values 
of ajl? . These plots clearly show that as a increases, the normalised gap frequency shifts to 
higher values. An additional feature is that the gap becomes softer with increasing a. These 
plots should be contrasted with figure O which shows similar conductivity plots, but now for 
scalar mass = — 3/Lg, and with both no backreaction and a backreaction of = 0.05. 
Three sample values of a are shown: 0, 0.125, and 0.25. In these plots, the gap if anything gets 
harder with increasing a, however once backreaction is included, the gap becomes more gentle 
and extended, and while the dip in iM(cr) is smoothed, it is still clearly apparent. 

In all of these plots, the gap is clearly indicated by a rise in the real part of a, which coincides 
with the global minimum of Im (cr). As already noted, the imaginary part of ()44p is only valid 
up to a linear term in u, the size of which is dependent on the renormalization scheme employed 
(and also on the charge Q). One can therefore tune this linear term in Im((t) to create a finite 
global minimum if it is not initially present, and indeed the plots in figure H] differ from those 




Figure 4. Conductivity: a range of plots showing the real (solid red line) and imaginary (blue 
dashed line) parts of the (normalized) conductivity as a function of (normalized) frequency for 
the scalar mass rn? = — 3/-L^. Plot (a) is a = 0, (b) is a = 0.1, and (c) is a = 0.2, and (d) is 
a = 0.25. 



presented in [Ij by precisely this sort of term so as to make the minimum apparent. This 
minimum can therefore be used to define ojg, the value of the frequency gap. 

The frequency gap is a distinct characteristic of a superconductor and in the BCS theory of 
superconductivity this frequency gap corresponds to the minimum energy required to break a 
Cooper pair. As mentioned above, in [11] it was claimed that for the holographic superconductor 
the relation LOgjTc — 8 had a certain universality, proving stable for a range of scalar masses 
and dimensions. In both [Ij and [2J, this relation was shown to be unstable to Gauss-Bonnet 
corrections. Figure [SKd) in particular gives a very clear indicator of how backreaction and 
higher curvature terms affect the gap. Increasing either a or increases ujg/Tc- For the case 
of increasing a, the effect occurs mainly because of a shift in the gap, rather than a significant 
alteration of Tc, which varies much more strongly with backreaction than a. On the other hand, 
varying practically does not alter ujg at all, whereas Tc drops dramatically, leading to a sharp 
rise in ujg/Tc- 

4. Summary 

This presentation reviewed the work of pLj and |2] on exploring the implications of Gauss-Bonnet 
corrections to holographic superconductors. The results show that increasing backreaction lowers 
the critical temperature of the superconductor hence increasing ujgjTc- The effect of higher 
curvature terms is more subtle. Although these initially act in a similar fashion to backreaction 
in lowering the critical temperature, for significant GB coupling and larger scalar masses, the 
critical temperature eventually begins to increase. The conductivity gap is also modified, with 
both iOg and altering to increase the ratio ujg/Tc- Clearly higher dimensional holographic 
superconductors have a rich structure, with or without higher curvature corrections. 
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Figure 5. Conductivity: a range of plots showing the conductivity of the = — 3/Lg 
superconductor. Plots (a) — (c) show the real (solid line) and imaginary (dashed line) parts 
of the conductivity as a function of frequency for a = 0,0.125, and 0.25 respectively. In each 
case the conductivity is shown for no backreaction in black, and for a backreacting parameter 
= 0.05 in magenta. The slight undulations in the plots at large a; is a numerical artefact. 
Plot (d) shows the gap frequency as a function of Tc with the line ujg = 8Tc is shown in black. 
The different colours represent from right to left: Red is = 0, Blue is = 0.05, Green 0.1, 
and Orange = 0.2. In each case a is incremented from to 0.25. As the gap alters rapidly 
near the Chern-Simons limit, the dotted lines are added by hand to guide the eye. 
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